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Abstract 



In previous papers we showed that in a one-armed deformed disks, p-mode and 

g-mode oscillations are resonantly excited by horizontal resonance, and applied it 

to high frequency QPOs observed in low mass X-ray binaries. In that model, the 

observed time variation of kHz QPOs is regarded as a result of a time-dependent 

q . precession of the deformation. In this paper we consider another possible cause of time 

5h ■ 

variation of kHz QPOs. That is, we demonstrate that in a two-armed deformed disks, 
p-mode and g-mode oscillations are excited by vertical resonance, not by horizontal 
resonance (horizontal resonance dampens them). Furthermore, we show that in the 
case of vertical resonance, the frequencies of disk oscillations excited can vary with 
■<^j- | time if vertical disk structure changes with time. A brief application of these results 

to the time variation of observed kHz QPOs is made. 

Key words: accretion, accrection disks — quasi-periodic oscillations — resonance 
^ . — neutron stars — two-armed disk deformation — X-rays; stars 
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1. Introduction 

Many authors now think that high frequency quasi-periodic oscillations (QPOs) ob- 
served in low mass X-ray binaries are some kind of disk oscillations in strong gravitational 
field, and that it gives a promising way to estimate the mass and spin of the central sources. 
High frequency QPOs both in black hole candidates and in neutron star sources will have the 
same dynamical origin (Abramowicz et al. 2003), although there are some differences in their 
frequency variations. 

Based on closeness to 3 : 2 of frequencies of twin QPOs, importance of resonant phe- 
nomena was emphasized by Abramowicz and Kluzniak (2001). Subsequently, from a different 
context, Kato (2004, 2008a, c) pointed out the importance of resonant processes in a deformed 
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(warped or eccentric) disks as an excitation process of high frequency QPOs. That is, non-linear 
interaction between disk oscillations and the deformed part of disks resonantly excite or dampen 
disk oscillations. There are two types of the resonance. One is a horizontal resonance (Lindblad 
resonance) and the other is a vertical resonance. The inertial-acoustic oscillations (p-mode) and 
the gravity oscillations (g-mode) are found to be excited by the horizontal resonance. 

The resonant excitation of the p-mode and g-mode oscillations by the horizontal reso- 
nance seems to occur most efficiently at the radius where the condition 

is satisfied (Kato 2004, 2008a,c), when the one-armed deformation of disks have no precession, 
where k(t) is the epicyclic frequency, Q(r) is the angular velocity of disk rotation and r is 
the radius from the central source on the disk plane. In the case of Schwarzschild metric, the 
resonance occurs at 4r g , where r g is the Schwarzschild radius defined by r g = 2GM/c 2 , M being 
the mass of the central source. The frequencies of the p-mode and g-mode oscillations excited 
by this resonance are 



where the subscript res denotes the values at the resonant radius derived from equation §1§ and 
m(= 1,2,...) are the wavenumber of the oscillations in the azimuthal direction. The frequencies 
specified by equation ([2]) have ratios of some rational numbers. 

The frequencies of oscillations given by equation (T5]) can qualitatively describe the high 
frequency QPOs in black-hole and neutron-star X-ray binaries. In the case of neutron-star 
X-ray binaries, however, there is an additional and challenging observational evidence that the 
frequencies of kHz QPOs change with time. If we want to describe this time change of kHz 
QPOs in the framework of the above-mentioned model, we must introduce the assumption that 
the one-armed disk deformation has time- dependent precession in the disks of neutron-star 
X-ray binaries^ , although such precession is not required in the black-hole accretion disks. If 
the frequency of precession of disk deformation is denoted by u p , the resonant condition where 
the horizontal resonance occurs is changed from equation (CQ) to 



The resonant radius, r res , changes with a time change of u pj and thus the frequencies of reso- 
nantly excited disk oscillations given by equation ([2]) changes with time with correlation. 

This may be one of causes of the time variation of kHz QPOs. However, a question 
naturally raised here is whether a variation of precession frequency that is enough to account 
for the time variation of kHz QPOs is really expected in disks of neutron-star X-ray binaries. 

1 In the case of neutron stars, the central sources have surfaces and the inner part of disks may have influences 
of stellar magnetosphere. In such situations, the warp might have precession (e.g., Meheut and Tagger 2009). 




(2) 
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Hence, it will be worthwhile considering other possibilities of frequency change of kHz QPOs. 
This is the purpose of this paper. 

First, we should notice that in the case of vertical resonance, possibility of frequency 
variation of resonant oscillations is high compared with in the case of horizontal resonance. This 
is because the radius where the vertical resonance occurs is sensitive to the vertical structure 
of diskj§. In spite of this interesting characteristics of the vertical resonance (Kato 2005), we 
already showed that the disk oscillations resulting from vertical resonance are damped, not 
excited, in the case where the disk deformation is one-armed (Kato 2004, 2008a, c). 

In the case where the disk deformation is two-armed, however, the situations are changed. 
That is, as is suggested in this paper, the resonance that can excite the p-mode and g-mode 
oscillations is not horizontal one, but vertical one. This means that if the disk deformation in 
neutron-star X-ray binaries is two-armed, different from the case of black-hole X-ray binaries, 
the observed frequency variations of kHz QPOs will be naturally accounted for by taking 
the viewpoint that kHz QPOs are p- and g-mode oscillations resonantly excited by vertical 
resonance. 

In section 2, we briefly summarize the essence of resonant excitation of disk oscillations in 
one-armed deformed disks as a preparation for examining the case where the disk deformation 
is two-armed. Resonant excitation of p-mode and g-mode oscillations by vertical resonance in 
vertically isothermal disks is examined in section 3. In section 4, the results in section 3 are 
extended to the case of vertically polytropic disks, and time variation of excited oscillations by 
a change of vertical disk structure is discussed. Numerical results are given in section 5, and 
the final section is devoted to discussion. 

2. Brief Summary of Resonant Excitation of Disk Oscillations in Warped Disks 

Before considering resonant excitation of disk oscillations in two-armed deformed-disks, 
we shall briefly summarize the essence of resonant excitation of disk oscillations in warped disks, 
since the procedure for examining the former problem is a simple modification of the latter. 

2.1. Unperturbed Disks 

For mathematical simplicity, the unperturbed disks are assumed to be isothermal in 
the vertical direction. This leads to the fact that the density is stratified exponentially in the 
vertical direction as 



where poo is the density on the equatorial plane (z = 0), and H(r) is the half-thickness of the 

2 In the horizontal resonance, the resonant radius is determined by the radial distributions of the Kcplcrian 
angular velocity and the radial epicyclic frequency, and they are almost time-independent in geometrically 
thin disks. 




(4) 
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disk, r being the distance from the rotating axis of the disks. 

The half-thickness, H, of disks is related to the vertical epicyclic frequency, Q±(r), by 



where c s is the isothermal acoustic speed. The vertical epicyclic frequency, Q±(r), is equal to 
the angular velocity of the Keplerian rotation, fl^(r), in the case of non-rotating central object, 
and practically equal to the angular velocity of disk rotation, Q(r). Hereafter, however, we do 
not use Qk or Q instead of Q± so that we can trace back the effects of Q± on the final results. 

2.2. One-armed Disk Deformation 

We assume that the disks described above are deformed from axisymmetric state by 
some external or internal cause. The deformation is a warp or an eccentric deformation in the 
equatorial plane. They are assumed, for simplicity, to be time-independent. 

The Lagrangian displacement associated with the deformation, £ w (r), is denoted by 



where if is the azimuthal direction of the cylindrical coordinates (r, ip, z), 7i n w is the Hermite 
polynomial of order n w with argument z/H. In the case of eccentric deformation in the 
equatorial plane, n w is zero (n w = 0), while it is unity (n w = 1) in the case of warp. It is noted 
that the number of nodes of £^ is smaller than those of ^ and ^ by one. This is true even 
for disk oscillations [see equation ((7j)]. 

2.3. Disk Oscillations 

Disk oscillations that are considered here are assumed to have moderately short ra- 
dial wavelength so that it is shorter than the characteristic length of radial variation of disk 
structure. Furthermore, the wave motions are assumed to occur isothermally. 

If the above approximations are adopted, the r- and z- dependences of disk oscillations 
on non-deformed disks are approximately separated (Okazaki et al. 1987). We can then express 
the displacement vector, £(r,t), associated with an oscillation mode of (u, m, n) as 



where u is the frequency of the oscillations, and 7i n is the Hermite polynomials as mentioned 
before. The integer n(= 0,1,2,...) specifies the number of nodes of £ r (and £^) in the vertical 
direction (^-direction). 

Hereafter, for simplicity, we characterize the oscillations and disk deformations by the 




(5) 



C W = exp(-z^)C w (r)^w(z/iJ), 
^ = exp(- i ^(r)HMz/H), 
C = eM-^(r)H nW ^(z/H), 



(6) 



£ r (r,t) =exp[i(tjt-m(p)]^ n (r)H n (z/H), 
^(r,t) =exp[i(wt — m<p)]£ Vf n{r)?i n (z/H), 
£ z (r,t) =exp{i(ut-mip)]£ Z:n (r)H n -i(z/H), 



(7) 
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set of (lj, m, n). For example, a warp belongs to a mode characterized by (0,1,1) and an 
eccentric deformation on the disk plane by (0,1,0). 

Here, classification of disk oscillation modes is briefly summarized (for details, see Kato 
2001; Kato et al. 2008). The oscillations with n = occur predominantly on the equatorial 
plane. They are inertial-acoustic oscillations and called hereafter "p-mode" oscillations. In the 
cases where n > 1, we have two kinds of oscillations. In one of them, (u — mQ) 2 < k 2 , while in 
the other one we have {uj — mQ) 2 > nQ\. We call the former "g-modes" , and the latter "vertical 
p-modes" except for some special cases mentioned below. In the case of n = 1 with m = 1 (and 
some special cases of n > 2 with m > 2), the latter oscillations are almost incompressible and 
have low frequencies. They are specially called "c-mode" oscillations. A warp is included to 
this type of oscillations. 

In the resonant excitation problem to be examined in this paper, we group the above 
various oscillation modes into two classes and treat the oscillation modes in each class as a pack, 
due to similarity of mathematical treatment. One class is p-mode and g-mode oscillations, and 
the other one is vertical p-mode and c-mode oscillations. 

2.4- Coupling between Deformation and Oscillations 

Nonlinear couplings between the deformation specified by equations ([6]) and the disk 
oscillations given by equations ([7]) induce disk oscillations that are characterized by (u, mil, 
h), where n = n i 1 when the disk deformation is a warp, while n = n when the deformation 
is an eccentric deformation on the equatorial plane. Arbitrary combinations of ± are allowed. 
We call the oscillations resulting from the coupling " intermediate oscillations" . 

The intermediate oscillations have resonant interaction with the disks at particular radii. 
One of resonances occurs at the radii where the intermediate oscillations with (u, mil, n) 
have the Lindblad resonances, which are specified by 



where k is the (horizontal) epicyclic frequency. This resonance is call hereafter "horizontal res- 
onance" . Another one occurs at the radii where the frequency uo of the intermediate oscillations 
of (a>, mil, ft) becomes equal to the eigen-frequency of vertical oscillations of disks. The radii 
are characterized by 



This resonance is callded hereafter "vertical resonance". 

2.5. Resonant Excitation of Oscillations 

The intermediate oscillations interact nonlinearly with the disk deformation, after having 
the resonance mentioned above, to feedback to the original oscillations. This feedback process 
amplifies or dampens the original oscillations. Detailed analyses (Kato 2004, 2008a, c) show 




(8) 




(9) 
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that in both cases of horizontal and vertical resonances, the growth rate of oscillations, —uj\ (oj\ 
being the imaginary part of frequency of oscillations), can be expressed in the form of 

sign[u;- (m± l)f2] rcs . . 

- oc , (10) 

where E is the wave energy of the original oscillations (uj, m, n) in consideration, and sign[c<j — 
(m±l)Q] res is the sign of uj — (m±l)fi at the resonant radius, r res jf| The value of the proportional 
coefficient of equation (fit)]) depends, of course, on modes of oscillations and types of resonances. 
However, in some typical cases, the coefficient is always negative definite (Kato 2004, 2008a, c). 

This means that the condition of excitation of disk oscillations is 
sign[o;-(m±l)^] rcs 

E 1 ' 

This condition allows us to have a simple physical interpretation. We notice first that 

for a resonance to occur efficiently, the resonant radius, r rcs , must be in the radial region where 
both the original and the intermediate oscillations are dominated. Furthermore, we notice 
that in general a wave with (ou, m) has a negative energy if the wave is dominated inside the 
corotation radius given by uj — mQ = 0, while it is positive if the wave is outside the corotation 
radius. This consideration suggests that we can regard the signfu; — (m± l)f2] res as the sign of 
the wave energy, E mt , of the intermediate oscillation. Furthermore, the sign of the wave energy 
of the original oscillation is the same as the sign of (uj — mf2) rcs . Based on these considerations, 
we can write the amplification condition (111)) as 

flint 



E 

or 

sign[u; — (m ± 



<0 (12) 
<0. (13) 



res 



sign(u; — mO) res 

This condition can be interpreted in the following way. If an oscillation with positive 
energy (E > 0) resonantly interacts, at a resonant radius, with an intermediate oscillation 
with negative energy (E mt < 0), both oscillations are amplified by energy flowing from the 
intermediate oscillation to the original oscillation. The original oscillation is amplified, since 
it has E > and receives positive energy. The intermediate oscillation also grows by loosing 
energy since E mt < 0. In the case of E < and E mt > 0, the resonance also amplifies the 
oscillations. In this case, the direction of energy flow is opposite: It flows from the original 
oscillation to the intermediate oscillation at the resonant radius. 

2. 6. Resonant Radius and Type of Resonance That Excites Oscillations 

It is noted that the resonant condition [equation (jSJ) or (Qj alone does not uniquely 
determine the radius of resonance. An additional restriction is necessary. We assume that the 



3 The resonant radius, r ros , is not uniquely determined by the resonant condition (|8|) or ([9]) alone. An 
additional condition is necessary, which will be discussed later in subsection 2.6. 
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original oscillations are most strongly excited when the resonance occurs near the boundary 
of their propagation region. Near the boundary the oscillations have long radial wavelength 
and their group velocity in the radial direction vanishes, i.e., they stay there for a long time 
compared with in other places, and thus will grow there most strongly. Hence, when we consider 
p- and g- mode oscillations, we assume that an additional condition to be adopted to specify 
the resonant radius is 

(uj-m{l) 2 -K 2 = 0, (14) 

since this represents the boundary of the propagation region of p- and g-mode oscillations. In 
the case where we consider vertical p-mode oscillations and c-mode oscillations, on the other 
hand, we assume that the resonance occurs at 

(u-mtt) 2 -ntt\ = 0, (15) 

by the same reason as the above. 

As mentioned in subsection 2.3, disk oscillations can be grouped into two classes in 
studying the present excitation problem. The first one is the p- and g-mode oscillations. The 
second class is the vertical p-mode and c-mode oscillations. Concerning the type of resonances, 
we have also two types, i.e., the horizontal and vertical resonances. Hence, we have four cases 
in combination of the set of oscillations and resonances. Among them, the excitation of disk 
oscillations occurs in the case where the oscillations are p- and g-modes and the resonance is 
horizontal (Kato 2004, 2008a, c). Combining equations (jHJ) and ( T141 . we find that the resonance 
in this growing case occurs at the radius where the condition of 

K = hn (16) 

is satisfied. The radius where this condition is satisfied is 4r g , i.e., r res = 4r g , when the metric 
is the Schwarzschild one. If the metric is the Kerr, the radius becomes smaller than 4r g . 



The frequencies of the disk oscillations that are excited there are [see equation (TT4"]) ] 

uj = (mfi ± ft) res . (17) 

Application of these resonantly-excited oscillations to high-frequency QPOs is made by Kato 
and Fukue (2006) and Kato (2008b). 

3. Resonant Excitation of Disk Oscillations in Two- Armed Disks 

After the above preparation, we now proceed to the main purpose of this paper, i.e., 
examination of the case where the disks are deformed from the axisymmetric steady state into 
a state with a two-armed pattern. Different from the case of a one-armed pattern, a two-armed 
pattern cannot be stationary in general, i.e., it will be time-dependent, and wavy. What we 
need here is that a pattern has approximately a constant frequency for a time interval longer 
than the characterisitic time by which disk oscillations grow. The origin of such a two-armed 
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pattern is a problem to be discussed and clarified, but here we simply assume that such a 
pattern exists on the disks by some internal or external causes. For example, numerical 3D 
MHD simulations of accretion disks (Machida and Matsumoto 2008) show that one-armed and 
two-armed patterns with slow rotation are produced in the innermost region of disks at a certain 
stage of disk evolution by magnetic-field streching and reconnection processes. 

The displacement vector, £ T (?",£), associated with a two-armed deformation is now ex- 
pressed as 

£(r,t) =ex V [ i (u; T t-2< f )}£j(r)H n T(z/H), 
ej(r,0 = exp[^ T t - 2ip)]^(r)HMz/H), 

g(r,t)=ex V ii(u T t-2 V )}i](r)H n T_ 1 (z/H). (18) 

Here, u T is the angular velocity of the pattern and taken to be a free parameter. As the integer 
n T , we are mainly interested in the cases of n T = 2 and 3 (see section 6). 

On such deformed disks, we impose disk oscillations. The displacement vector, £(r,t), 
associated with the oscillations is described again by equation ([7]). 

3.1. Resonant Conditions 

The nonlinear coupling between the disk deformation described by (u; T , 2, n T ) [see 
equation ( fi~8l) ] and the disk oscillations described by (u, m, n) [see equation (0)] induces the 
intermediate oscillations of (u;±c<; T , m±2, n), where n = n±n T . These intermediate oscillations 
have resonances with the disk rotation. The horizontal resonance occurs at [cf., equation ([8])] 

[u±u' T -(m±2)Q] 2 -k 2 = 0, (19) 

and the vertical resonance occurs at [cf., equation ((9])] 

[cu±cu T -(m±2)fi] 2 -n^ = 0. (20) 

As mentioned in subsection 2.6, an additional condition is necessary to determine 
uniquely the radius of resonance. We focus our attention again to the case where the reso- 
nance occurs at the radius where the group velocity of the original oscillations with (a;, m, n) 
vanished to stay there for a long time. That is, when we consider p- and g-mode oscillations, 
we adopt [equation (fT4|) ] 

(u)-mtt) 2 -K 2 = 0, (21) 

as the additional condition. In the case where the excitation of vertical p-mode and c-mode 
oscillations is examined, we adopt [equation ( |T5l) ] 

(uj - mn) 2 - nn\ = 0. (22) 

There are four cases in combination of types of resonance and types of oscillation, i.e., 
two cases (horizontal or vertical resonance) for p- and g-mode oscillations, and two cases (hor- 
izontal or vertical resonance) for vertical p-mode and c-mode oscillations. Before examining 
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these cases separately in subsection 3.3, we consider the excitation condition. 

3.2. Excitation Condition 

In the case of one-armed deformation of disks, the condition of excitation of disk os- 
cillations is given by equation ([TO]) . The mathematical procedures to derive the condition are 
complicated (see Kato 2008a, c), but they can be straightly extended to the case of two-armed 
deformation. The results show that the growth rate of oscillations, — U[, can be expressed in 
the form of 

sign[u; ± u p - (m ± 2)Q] res 

UJ\ OC — . 

That is, compared with the case of the one-armed deformation, in the present case of two-armed 
deformation, uj is changed to ui ± u p and m ± 1 is changed to m ± 2. 

The next problem is to examine the sign of the proportional coefficient on the right 
hand side of equation ( 1231) . A straightforward generalization of the procedures of the one- 
armed deformation suggests that in some simplified cases (for example, the case where the non- 
linear coupling terms between the original oscillation and the deformation are constant in the 
resonant region^), the proportional coefficient is negative definite as in the case of one-armed 
deformation. When the proportional coefficient is negative definite, we have again a simple 
physical interpretation of equation ( 1231) [see subsection 2.5]. Considering them we suppose that 
the amplification condition in the present two-armed case is given by 
signfcj ± uj p — (m ± 2)0] r 



J res 



< 0. (24) 
signfw — rmzjres 

In other words, the condition is also expressed as E mt / E < 0, and has a simple physical meaning. 

3.3. Growing Cases and Their Resonant Radius 

Based on the resonant conditions [relevant combinations of equation (TTTJ|) or (12"U1) to 
equation fT2TT) or fT22l) ] and the excitation condition [equation (T2"4"|) ]. we now examine what types 
of oscillations (p- and g-mode oscillations or vertical p- and c-mode oscillations) are excited by 
what type of resonances (horizontal or vertical resonance) and where the radii are. 

i) Horizontal resonance of p- or g-mode oscillations 

In this CBjSC ; clS mentioned in subsection 3.1, the resonant radii are the places where both 
equations ffl9|) and (T2Tj) are simultaneously satisfied. More explicitely, horizontal resonances of 
p- and g-mode oscillations occur at the radii where one of the following set of two equations 
are satisfied: 

(a): u ±lu t — (m±2)il = k and uj — mQ, = k, 



The resonant region has a finite width in the radial direction around the resonant radius. The width 
depends on the disk temperature. In a pressureless disks, the width is infinitcsimally narrow, but increases 
with increase of disk temperature (e.g., Kato 2008c). 
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(b) : uj ± uj t — (m ± 2)tt = k and uj — mVt = — k, 

(c) : uj ±uj t — (m±2)£l = —k and uj — m£l = K, 

(d) : uj ± uj t — (m ± 2)fi = — k and co> — mf2 = — k. (25) 

By inspection we see that the radii satisfying condition (a) or condition (d) are not 
interesting here, since the excitation condition ( 1241) is not satisfied in these cases. If there 
are radii where condition (b) or (c) is satisfied, the oscillations are excited there since the 
excitation condition (124"1) is satisfied there. From condition (b) or (c) we see that such radii 
are ±uj t =p 2Q = 2k. Since uj t will be much smaller than Q in practical cases, we take as the 
resonant radius where oscillations are excited 

« = n±^-, (26) 

where + is for the case of uj t < and — is for uj t > 0, since k < Q. In the case where uj t is 
much smaller than Q, this resonance occurs at an outer region of the disks, and the frequency 
ratios of excited oscillations are roughly 1 : 2 : 3.... Such oscillations will be of interest, but 
are subjects outside the present issue. 

ii) Vertical resonance of p- or g-mode oscillations 

In this case, the radii where resonance occurs efficiently are places where both equations 
(1201) and ( T2TT) are simultaneously satisfied. These conditions can be written down in the following 
four cases: 

(a) : uj ± uj t — (m ± 2)Q = h 1 ^ 2 ^^ and uj — mVL = n, 

(b) : uj ± uj t — (m ± 2)f2 = n 1//2 f2j_ and uj — mVt = — k, 

(c) : uj ± uj t — (m ± 2)Q = — n 1//2 f2j_ and uj — mVt = K, 

(d) : uj±uj t -(m±2)n = -n 1/2 n ± and uj-mVt = -K. (27) 

As in the horizontal resonance, we see that cases (a) and (d) are uninteresting even if 
they have solutions, since the excitation condition ff2^|) is not satisfied. Resonances resulting 
from case (b) or case (c), on the other hand, satisfy condition (12^1) . In cases of (b) and (c), the 
resonant radii are found to be 

K = 2tt-n 1/2 {l ± Tu T , (28) 

where both signs of ± are possible. The negative sign is for the case of (b) and positive one 
is for (c). Here, uj t has been assumed to be much smaller than Q. This resonant condition is 
satisfied in the case of n = 2 and 3 at inner region of relativistic disks. 

If the resonant radius, r res , is determined by solving equation (j28|) . the frequencies of 
resonantly excited oscillations are found to be 

UJ = {mVL ± ft) res - (29) 
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iii) Horizontal resonance of vertical p-mode or c-mode oscillations 

In this case, the set of equations to be used to determine resonant radii are [see equations 
dm and Ol 



(a) : u ±u T — (m±2)Q = n and u — mVt = n 1 ^ 2 fl±, 

(b) : lo ± uj t — (m ± 2)Q = k and lo — mVt = — ra 1//2 Oj_, 

(c) : u ± u T — (m ± 2)Q = — hi and lo — mVt = n 1 ^ 2 Q,±, 

(d) : uj±u t - (m±2)Q = -k and u - mVt = -n 1/2 £l ± . (30) 

As in the previous studies of i) and ii), cases (a) and (d) are outside of our present interest, since 
the oscillations that satisfy the conditions do not satisfy the excitation condition (}24l) . In the 
cases of (b) and (c), on the other hand, oscillations are excited, and the radius is characterized 
by 

k = 2Q - n 1/2 Q ± ± u T . (31) 

This expression for resonant radii is the same as relation (j28p . except that n in relation (I28p is 
now replaced by n. 

The frequencies of resonantly excited oscillations are 

uj = (mtt±n 1/2 tt ± ) ms , (32) 

where the resonant radius, r res , is now determined by equation ( 13T1) . 

Finally, it is noted that there is no resonance characterized by the set of equaions (TSUI) 
and ( l22i) . as far as n and n are moderate integers. 

4. Frequency Variation by Change of Vertical Disk Structure 

Among resonantly excited oscillations discussed in the previous section, the p- and g- 
mode oscillations resulting from vertical resonance [i.e., the set described by equations (T28l) 
and ( 1291) ] are of interest in relation to the observed kHz QPOs, since their frequencies are in a 
reasonable frequency range, as discussed in previous papers. 

If we want to account for the observed frequency variation of kHz QPOs by the model 
described by equations (1281) and (1291) . the frequency of two-armed pattern, lj t , must change with 
time. Variation of u T changes the resonant radius [see equation (128!) ] and thus the frequencies 
of oscillations [see equation ( 1291) ] are changed. This may be one of possible causes of frequency 
changes of the observed kHz QPOs, but it is not clear whether a large variation of u T required 
to explain the time variation of kHz QPOs is generally expected. 

Here, one of another possibilities of time variation of resonant radius is considered. So far, 
we have assumed that the disk is vertically isothermal and oscillations also occur isothermally. 
Now we relax this assumption, and consider the case where the pressure, p, and density, p, 
are distributed in the vertical direction with a polytropi relation, i.e., p = Kp l+1 I N , and the 
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polytropic index N changes with time. It is noted that in the polytropic disks the vertical 
integration of the vertical hydrostatic balance gives 



z 2 



T (r,z)=Too(r)(l-— J, (33) 



Z 2^N 



Po{r,z)= Pm {r)[l-—j , (34) 

p Q (r,z)=T 00 (r)(l-—) , (35) 

where subscript represents the quantities in the equilibrium state and 00 are those on the 
equatorial plane (e.g., Kato et al. 2008). 

We consider adiabatic disk oscillations on such polytropic disks, assuming that the ratio 
of the specific heat, 7, is related to N by 7 = 1 + 1/N. Compared with the case of isothermal 
disks, derivation of dispersion relation of oscillations in such polytropic disks are much compli- 
cated (e.g., Perez et al. 1997; Silbergleit et al. 2001). However, in the limiting case where the 
couplings between horizontal and vertical motions in an oscillation mode are neglected, we can 
easily see that the eigen-frequency of local vertical oscillations with (u, m, n) is given by 

(u-mtt) 2 -^ n tt 2 ± = 0, (36) 

where for the fundamental (n = 1), the first overtone (n = 2) and the second overtone (n = 3), 
we have (e.g., Kato 2005) 

*i = l, 

^2 = 2 + ^ = 1 + 7, 

*3 = 3 + ^ = 3 7 . (37) 

It is noted that in the case of isothermal disks, 1/N — and \l/ n = n. 

The above consideration suggests that the condition of vertical resonance, equation (120]) . 
is now changed to 

[u±u' T -(m±2)Q} 2 -ty n Q 2 ± = 0. (38) 

This relation shows that the resonant radius of the vertical resonance depends not only on u T 
but also on 7. In real accretion disks, a change of mass accretion rate, for example, may bring 
about a change of the disk vertical structure. This change gives rise to a change of radius of 
the vertical resonance, leading to frequency change of resonantly excited disk oscillations. 

The next subject to be investigated is how the excitation condition (j24l) is changed 
in the present case of polytropic disks. Detailed investigation of this is very complicated in 
mathematical treatment, since normal mode analyses of oscillations are troublesome in the 
case of polytropic disks. Let us first consider the separability of variables associated with 
oscillations. If the radial variation of oscillations is strong, a physical quantity associated 
with the oscillations, say V(r, z), is approximately separated as V r {r)V z (z). In the case of 
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isothermal disks, V z (z) is the Hermite polynomials [see equation 0], while it is the Gezenbauer 
polynomials in the case of polytropic disks [see Perez et al. 1997; Silbergleit et al. 2001]. This 
makes description of oscillations complicated. 

Furthermore, when we want to describe the non-linear couplings between oscillations, 
we must expand a product of Gezenbauer polynomials into a series of the Gezenbauer polyno- 
mials, and must use orthogonal relations among the polynomials to separate oscillation modes. 
These are much complicated compared with in the case of the Hermite polynomials. These 
make a detailed examination of excitation condition troublesome, compared with in the case of 
isothermal disks. Hence, a detailed derivation of excitation condition in the case of polytropic 
disks is beyond the scope of this paper. Here, we must be satisfied with a rough and physical 
considerations. That is, we assume that the excitation condition ( J24l) still hold even in the case 
of polytropic disks, since it has a simple and reasonable physical meaning that will be free from 
a particular disk structure, as mentioned before. 

Hereafter we focus our attention only on the case where p- and g-mode oscillations are 
excited by the vertical resonance. In this case the basic equations to be used to specify the 
resonant radius are the set of equations (l38l) and (l29il . and the equation to be used to judge 
whether the oscillations are really excited is equation (1241) . 

By generalizing the procedures in ii) in subsection 3.3, we easily see that the resonant 
radii, r res , in the cases where oscillations are really excited are obtained by solving 

K = 2n-^l /2 n ± T^ T , (39) 

and the frequencies of the oscillations, u>, are given by 

UJ = {mVt ± ft) r es- (40) 

5. Numerical Results 

As is shown in equation ( l40l) . the frequencies of oscillations excited at the resonant radius 
are a discrete set, characterized by m and ±. As the typical frequencies, we take here (Q — ft) res 
and (2ft the previous papers. The former oscillation has m — 1, while the latter 

does m — 2. We take a picture that the observed QPOs come from high energy photons that 
are Comptonized in a corona surrounding a geometrically thin disk where the oscillations are 
generated. If this picture is adopted, the oscillations with m — 1 are observed in the twohold 
frequency (Kato and Fukue 2006). Based on this situation, we consider that 2(Q — k) vcs and 
(2Q — K)res correspond to the typical twin frequencies of the observed kHz QPOs, and denote 
them as 

2u LL = 2(Q - K) rcs , u L = (2Q - n) res . (41) 

The figures of frequencies 2u> LL and u>l depend on the resonant radius, r res , but the 
relation between 2u;ll and u;l is free from detailed models determining r rcs . That is, the relation 
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Fig. 1. Diagram showing the relation between cjl and 2cjll for some values of mass, M, and spin pa- 
rameter, a*. The set of M (in units of Mq) and a* adopted to draw the curves are, from the uppermost 
to the lowermost curves, (2.0, 0.2)(green), (2.0, 0.0) (red), (2.4, 0.2) (blue), and (2.4, 0.0)(violct). Diagram 
showing the frequency correlation of the twin kHz QPOs of some neutron-star X-ray sources (taken from 
Abramowicz 2005) arc superposed by assuming that the upper and lower kHz QPOs correspond, respec- 
tively, to wl and 2wll- 

depends only on the mass of the central object, M, and the spin parameter, a*, representing 
the metric. The variations of such parameters as 7 and u T determine only the allowed range 
of variation on the c^l - 2u;ll curve. The u;l - 2cl>ll relation is shown in figure 1 for some sets of 
M and a*. On this figure, the diagram showing the observed relation between the upper and 
lower frequencies of the twin kHz QPOs of some typical sources has been superposed, assuming 
that the lower and upper kHz QPOs correspond, respectively, to 2cl>ll and wl- In the region of 
2wl L > 600 Hz, the - 2o>ll relation seems to well describe observations if M = 2.4M and 
a* = are adopted. 

Next, we examine which part of the u;l - 2u;ll curve is allowed in the case where 7 and 
up- vary in resonable ranges. As a preparation to this study, we first examine how the resonant 
radius varies by changes of 7 and u . We consider first the case of u T = and h — 2. From 
equation (1391) . we see that in this case the resonant radius is described by 



Figure 2 shows how the resonant radius described by equation (|42j) changes as a function of 7 
for two cases of a* = and a* = 0.3. The change of resonant radius by change of u T is shown 
in figure 3 by adopting the positive sign in equation ( j39l ): 



In figure 3, 7 = 4/3 has been adopted and two cases of a* = and a* = 0.3 are shown. The r 




(42) 




(43) 



- u T relation in the case of k = 2Q — (1 + 7) 
in figure 3. 



1/2 _ 



u T is obtained by just changing the sign of o; T 
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Fig. 2. r res /r g - 7 relation obtained by solving equation (|42|) . w T = has been adopted. Two cases of 
a* = and a* = 0.3 are shown. 
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Fig. 3. r res /r g - w T relation obtained by solving equation ([4*3")) . Two cases of a* = and a* = 0.3 are 
shown with 7 = 4/3. 

We next examine how much 2cjll and cul vary, when u T and 7 are changed in reasonable 
ranges. Two cases of 7 = 1 and 7 = 4/3 are considered, and for each case of 7, u> T is changed 
from to -80 Hz (retrograde precession of deformation), with M = 2.4 M & and a* = 0.0. The 
results are shown in figure 4, where and l} t are shown as functions of 2co>ll- in the Z-sources, 
the frequency of the horizontal branch oscillation (HBO) is known to change with correlation 
with the frequencies of kHz QPOs, i.e., the frequency of HBO is roughly (1/15) of the lower 
kHz QPO frequency. So, in figure 4 a straight curve of (1/15) x 2u;ll has been added in order 
to compare it with the |u; T | - 2u;ll relation obtained here. 

The comparison of the line of (1/15) x 2o>ll with the |u; T | - 2u;ll curve in figure 4 suggests 
that if changes of a> T and 7 are not independent but correlated, the present model may describe 
the observed frequency correlation between kHz QPOs and the horizontal branch QPOs (see 
the next section). 
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Fig. 4. Diagram showing the wl - 2wll and lo 1 - 2cjll relations in the cases where 7 = 1 and 4/3, in 
each case w T being changed in the range of Hz to - 80 Hz (retrograde precession). For a comparison, the 
curve of (1/15) x 2ojll is shown in order to compare with observations, since the frequency of horizontal 
branch oscillations and that of the lower kHz QPOs are correlated with the former being about 1/15 of 
the latter. 

6. Discussion 

Comparison of the present disk oscillation model of QPOs with the observed kHz QPOs 
suggests that masses of the central sources are around 2.4M Q , if they have no spin (see figure 
1). If they have spin, a larger mass is required to describe observations by our model. This 
mass required by our model seems to be rather large compared with that usually supposed as 
neutron star mass, although it is not excluded theoretically. 

So far as we assume that 2(£1 — /t) rcs (= 2c<jll) and (2a; — n) TCS (= correspond, respec- 
tively, to the upper and lower kHz QPOs, the above conclusion seems to be robust, since the 
2o;ll - &l relation does not depend on detailed models of resonance. This is true even in the 
case of one-armed deformed- disks, as far as 2u; LL and u LL are assumed to correspond to the 
observed twin kHz QPOs. 

There are a few possibilities to evade the above conclusion, since some assumptions and 
simplifications are involved in the present model. First, we have assumed that the angular 
velocity of disk rotation, Q, is the Keplerian. In the central part of the disks, this simplification 
might be violated, since magnetic field anchored to the central source may be strong enough to 
modify the disk rotation. 

More importantly, we should remember that our analyses are based on an idea that 
QPOs are propagating transient phenomena. That is, in addition to the condition of vertical 
resonance, equation (1381) . we imposed a condition that the oscillations that are excited most 
strongly are those whose radial group velocity just vanishes at the resonant radius, i.e., equation 
( 12T1) . Combining these two relations, we have obtained the resonant radius, equation ( 1391) . The 
imosed condition ( l2~TT) . however, is not always clear whether it is most relevant to specify the 
resonantly excited oscillations. If this condition is relaxed, the resonant radius given by equation 
( 1391 ) is modified. For example, the resonant radius may differ for each oscillation mode, and 
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2u>ll - relation is changed. Furthermore, we should notice that if the kHz QPOs are trapped 
oscillations, their frequencies are determined by a trapping condition, rather than by vanishing 
of group velocity. Then, a different approach from the present paper is necessary to consider 
frequencies of QPOs (see Ferreira and Ogilvie 2008; Oktariani et al. 2009). 

In this paper we considered the case of n = 2, i.e., the vertical component of displacement 
vector associated with the intermediate oscillations has two nodes in the vertical direction. In 
the case of n = there is no vertical resonance, and in the case of n — 1, the frequency of 
the vertical oscillation and that of the vertical resonance are always Q±, independent of the 
vertical structure of disks. Hence, n = 2 is the possible smallest value of n where the vertical 
disk structure can affect the frequency of resonant oscillations. The intermediate oscillations 
with n = 2 are realized when the set of (n, n T ) is (0,2), (1,1), (1,3), or (2,0),... in the case of 
vertically isothermal disks. We assume that the situation does not change much even in the 
case of polytropic disks. Then, if we remember that we are now treating the deformation with 
m T = 2, the sets with n T = 2 or n T = 3 are of interest among the above sets of (n, n T ), since in 
the cases of (m T = 2, n T = 2) and (m T = 2, n T = 3), we can expect low frequency deformation 
of disks by the following reason. 

Silbergleit et al. (2001) examined global adiabatic disk oscillations in polytropic disks, 
corresponding general considerations of the process deriving relations ( 1371) . By separating 
approximately a disturbance associated with disk oscillation, say V(r,z), into a separated form, 
say V r (r)V z (z) , they solved the resulting wave equation by WKB methods. They show that the 
frequency of vertical oscillation, say uj, is given by 

(uj-mtt) 2 = -n[j(n - 1) + (3 - n)]nl, (44) 

where 7 is the ratio of specific heats, m is the azimuthal wavenumber and n = 0, 1,2, ... is an 
integer characterizing the node number in the vertical direction. In the cases of (m = 2, n = 2) 
and (m = 2, n = 3), the above equation has solutions of the forms [see equation (137j) ] : 

u = -( 1 + l) 1/2 n ± + 2n in = 2) 

u = -(3 7 ) 1/2 fi ± + 2fi (n = 3). (45) 

Here, the above results are applied to disk deformation, and thus u and n in equation 
(|45l) are regarded, respectively, as u T and n T . Equation fj45|) then suggests that a low frequency 
disk deformation is expected when n T = 2 and n T = 3 for relevant figures of 7. For example, 
in the case of n T = 3, a steady deformation of disks is possible for 7 = 4/3, and u T is negative 
(retrograde precession) for 7 > 4/3, |u; T | increasing with increase of 7. This 7 - dependence of 
|o; T | is qualitatively the same as the u T - 7 relation required to describe observational correlation 
between kHz QPOs and HBOs. For example, let us consider the case where u T is related to 
7 as uj t = —100(7 — 0.5) (retrograde precession), and 7 changes in the range of 7 = 2/3 to 
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Fig. 5. Diagram showing 2wll - and \oj t \ - ojl relations in the case where uj t and 7 is related as 
lu t = —100(7 — 0-5) and 7 varies in the range of 2/3 to 4/3. The mass of the central source and the 
spin parameters are, respectively, M = 2.4M Q and a* = 0. The diagram has been superposed by the 
diagram showing the observed frequency correlations among kHz QPOs and low frequency QPOs (taken 
from Boutloukos et al. 2006), assuming that cjl corresponds to the higher kHz QPO. 

7 = 4/3i The 2uj ll - uj l and |u; T | - c^l relations in this case are shown in figure 5. On this 
figure, observed QPOs data for typical sources are superposed, assuming that cul corresponds 
to the upper kHz QPOs. 

Finally, we should emphasize that in our QPO model, non-linear couplings and reso- 
nances that are considered are between disk deformation and oscillations. Concerning oscil- 
lations themselves, however, our resonant model is linear; non-linear and resonant processes 
among oscillations themselves are not considered. If we want to describe the fact that the 
observed amplitudes of neutron-star twin QPOs changes sign as the observed frequency ratio of 
the QPOs passes through the value 3 : 2 (Torok 2009), non- linear resonant processes between 
twin QPOs should be considered as Horak et al. (2009) did and succeeded to describe it. The 
non-linear resonant processes among oscillations, however, are not considered in our resonant 
model, since in our model they are not main processes for determining the oscillations excited 
and their frequencies. 

5 We do not insist that uj t and 7 should be correlated in this way. This is just an example. 
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